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Abstract: A neutralino relic density in accord with WMAP measurements can be found 
in the minimal supergravity (mSUGRA) model in several regions of parameter space: the 
stau co-annihilation corridor at low m,Q and the hyperbolic branch/focus point (HB/FP) 
region at large mo at the edges of parameter space, and the bulk and Higgs boson funnel 
regions within. In the regions at the edge, the mass gap between the next-to-lightest SUSY 
particle (NLSP) and the lightest SUSY particle (LSP) becomes small, and backgrounds 
from 77 — ► ff (/ is a SM fermion) become important for NLSP detection at an e + e~ 
linear collider. We evaluate these backgrounds from bremsstrahlung and beamstrahlung 
photons, and demonstrate that these do not preclude the observability of the signal for the 
cases where either the stau or the chargino is the NLSP. We also delineate the additional 
portion of the stau coannihilation region, beyond what can be accessed via a search for 
selectrons and smuons, that can be probed by a search for di-tau-jet events plus missing 
energy. The reach of a LC for SUSY in the HB/FP region is shown for an updated value 
of rrit — 180 GeV as recently measured by the D0 experiment. 
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1. Introduction 



Weak scale super symmetry is a highly motivated extension of the Standard Model (SM)(T|. 
In models with gravity-mediated supersymmetry breaking and .R-parity conservation, the 
lightest SUSY particle is usually the lightest neutralino Z\, which is absolutely stable and 
serves as a good candidate particle for cold dark matter (CDM) in the universe^]. Indeed, 
the recent precision mapping of anisotropics in the cosmic microwave background radiation 
by the WMAP collaboration^, ||] has allowed a determination of 

VcDAih 2 = 0.1126±g;8iii (1.1) 

at the 2a level, where f2 = pcdm/Pc is the ratio of the cold dark matter density to the 
critical mass density of the universe, and h is the scaled Hubble constant. This DM is 
unlikely to all be massive black holes (at least in our galactic halo) ||, or ordinary baryons 
(in the form of brown dwarfs) since a baryonic density of this magnitude would lead to 
conflicts with both Big Bang nucleosynthesis and the density of baryons as determined 
from the acoustic peaks in the CMB spectrum. 

The minimal supergravity model (mSUGRA) forms a convenient template for explor- 
ing the experimental consequences of gravity- mediated SUSY breaking models Q. The 
mSUGRA model is characterized by four SUSY parameters together with a sign choice, 

mo, fn\jii Aq, tan/? and sign{p). (1-2) 

Here tjiq is the common mass parameter for all scalar particles at Q = Mqut, m i/2 is the 
common gaugino mass at Mqut, A) is the common trilinear soft term at Mqut, tan/3 is the 
ratio of Higgs field vacuum expectation values at the scale Mz, and finally the magnitude 
- but not the sign - of the superpotential [i term is determined by the requirement of 
radiative electroweak symmetry breaking (REWSB). In addition, the top quark mass rrit 
must be specified. While the PDG quotes a value nit = 174.3 ± 5.1 GeV[f?J, we note that 
the D0 experiment has recently determined a value mt = 180.1 ±5.1 GeV by a re-analysis 
of Run 1 data[§|. 

Over the years, several groups have evaluated the relic density of neutralinos in the 
context of the mSUGRA model|| [Hi. This has been recently re-examined in light of the 
WMAP data|ll]]. A value of ^ Zl h 2 in accord with the WMAP determination can be found 
in four broad regions of model parameter space: 1 

• The bulk annihilation region occurs at low mo and low mi/2 for all tan/3 values. 
In this region, neutralino annihilation in the early universe occurs dominantly via 
t-channel exchange of light sleptons. 

• The stau co-annihilation region occurs at low niQ where ~ m z ■> an d where 

fx — Z\ and f\ — f\ annihilation in the early universe serve to reduce the neutralino 

relic density to sufficiently low values |lO|l . 

1 We note that the CDM may consist of several components, so that WMAP observation should be 
interpreted as an upper bound f^fo 2 < 0.129 on the neutralino relic density. If we make the additional 
assumption that relic neutralinos consititute all the CDM, more stringent constraints are possible. 
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the hyperbolic branch/focus point region 

U, U, (HB/FP) at large mo near the 
edge of parameter space where \x becomes small, and the Z\ has a significant higgsino 



component which facilitates a large annihilation rate[15, 16, 17]. The location of this 
region is very sensitive to the value of mt pig ] . 

• the A-annihilation funnel at large tan/3 where mu and rriA — 2m % and ZiZi — > 
A, H —> ff (/'s are SM fermions) through the very broad ^4 and -ff resonances]^. 

The bulk annihilation region is now disfavored and possibly excluded because it gives rise 
to values of in violation of LEP2 limits, and can give large deviations from BF(b — ► 57) 
and (5 — 2) tt fi"cfl . A portion of the bulk annihilation region may survive, but only where it 
overlaps with the stau co-annihilation region or the A annihilation funnel. Other regions of 
parameter space may also give a reasonable relic density, including the stop co-annihilation 
region (for very particular Aq values) p0|] and the h annihilation corridor, where 2m^ — frih- 
Once the dark matter allowed regions of mSUGRA parameter space are identified, it 
is useful to see what the implications of these regions are for collider experiments. The 
reach of the Fermilab Tevatron]^, |2l| and CERN LHC|2^] have recently been worked out 
for the DM allowed regions of parameter space. In addition, the reach of a linear e + e~ 



collider has also been investigated for the DM allowed regions of the mSUGRA model) 23]. 



In Ref. [23], it was found that a yfs = 500 — 1000 GeV LC could probe the entire stau 



co-annihilation region for tan (3 ~ 10 by searching for selectron and smuon pair events. 
For higher tan/3 values, the stau co-annihilation region increases in mirt, while the reach 
for dilepton pairs decreases. The decrease in dilepton reach occurs mainly because the 
large-reach region is subsumed by the unallowed region where the f\ becomes the LSP (a 
fi LSP is not allowed in ii-conserving models because it would lead to charged stable relics 
from the Big Bang, for which there are stringent search limits). Some additional reach may 
be gained in the stau co- annihilation region by searching for stau pair production events, 
although this possibility was not investigated in Ref. |E3] . 



The reach of a LC for SUSY in the HB/FP region was also examined in Ref. 23]. It 
was found that in this region a chargino pair production signal could be seen above SM 
background essentially up to the kinematic limit for production of two charginos at the 
LC. Standard cuts for identifying I + 2 — jet events allowed chargino detection in the low 
m i/2 portion of the HB/FP region. In the high m^j portion of the HB/FP region, the 
W\ — Z\ mas gap became so small that specialized cuts were needed to pick out the signal 
from the background. In addition to usual SM backgrounds from the production of heavy 
particles (W, Z, t), contributions to backgrounds from e + e~ — > ff and 77 —►//(/ = b, c) 
were also included. The two photon cross sections were evaluated using Pythia [p4| for 
the bremsstrahlung photon distribution, while the beamstrahlung contribution was not 
included. It was found that in the HB/FP region, an e + e~ LC with yfs = 500 — 1000 GeV 
could have a reach for SUSY considerably beyond that of the LHC! 

It is significant that a neutralino relic density in accord with WMAP analyses occurs 
along the boundary of parameter space for both the stau co-annihilation case and the 
HB/FP region. In both these cases, the next-to-lightest SUSY particle (NLSP) becomes 
nearly degenerate in mass with the LSP. Thus, although NLSPs may be produced at 
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large rates at colliders, NLSP — > LSP decays may result in only a small visible energy 
release, potentially making their detection difficult. This is because 77 — > ff production, 
where the photons may originate from bremsstrahlung off the initial state electrons, or 
via beamstrahlung constitute an important background to the signal. In the stau co- 
annihilation region, we do not expect the two photon background to be problematic for the 
search for selectron or smuon pairs. The reason is that the e + e~ or ^ + /U~ pair originating 
from 77 annihilation is expected to come out back-to-back in the transverse plane for 
the background, but not so for the SUSY signal. However, for stau pair production, the 
77 — > tT background can be important, since the tau lepton is unstable, and the visible 
tau decay products will in general not be back-to-back in the transverse plane, especially 
when the daughter taus are relatively soft. 

One of the purposes of this paper is a careful examination of the two photon back- 



ground, including both bremsstrahlung and beamstrahlung contributions [25]. We describe 
our treatment of the 77 — ► // background , and its inclusion into Isajet v7.70[26] in Sec. |2[ 
We note here that it is likely that far forward detector elements reaching as close as 25 
mrad of the beam line will likely be implemented in a realistic detector. The main purpose 
of the forward detectors will be to veto 77 or events where the initial e suffers a slight 
deflection into the instrumented regions. In this work, we operate under the assumption 
that the initial state e are collinear with the beam. In this sense, our results may be 
regarded as conservative. 

In Sec. |3[ we examine extending the LC reach for SUSY in the stau co-annihilation 
region by searching for signals with two tau jets plus missing energy [ 27 1. The additional 



reach gained beyond that from the di-electron or di-muon channel is in fact substantial, 
especially for large tan/?. In Sec. ||, we re-examine the reach for SUSY by a LC in the 
HB/FP region, including the 77 — > bb and cc background to £ + 2-jet events originating 
from both bremsstrahlung and beamstrahlung processes. The added background from 
beamstrahlung is substantial, but can be eliminated by invoking an additional angular cut 
in addition to cuts already suggested in Ref . |23|] . We present updated reach results for one 
value of tan (3 at a higher value of mt = 180 GeV as suggested by a recent D0 analysis |§. 
In Sec. H we present a summary and conclusions of our results. 

2. Two photon backgrounds 

Our goal in this section is to evaluate the e + e~ — ► // background including contributions 
from bremsstrahlung and beamstrahlung, and including two photon annihilation processes. 
In this case, the e + e~ — ► // cross section is given by 

a(e + e~ ->■ //) = J dx x dx 2 [f e /e(xi,Q 2 )fe/e(x2,Q 2 )v(e + e~ -> //) 
+ Q 2 ) + /5T(*l)) {fj/T(^, Q 2 ) + ffiTfaj) a( 77 - //")] • (2.1) 

Here, f e / e (x,Q 2 ) is the parton distribution function (PDF) for finding an electron inside 
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the electron, and is given by the convolution 

■i 



fe/e(x,Q 2 )= dzf^ m (-,Q 2 )f b r m (z)/z, (2.2) 



x 



where f brem (x, Q 2 ) is the bremsstrahlung parton distribution function of the electron, and 
f beam (x) is the beamstrahlung parton distribution function of the electron. We use the 
Fadin-Kuraev distribution function for bremstrahlung electrons [28], given by 



f h e )t m (x, Q 2 ) = \P(1- x)^ 1 (l + |/9J - \P (1 + x) , (2.3) 

where f3 = ^(log^j — 1). For the beamstrahlung distribution function of electrons in 
the electron, we use the parametrization of Chenp9|, which is implemented in terms of 
the beamstrahlung parameter T (which depends in various characteristics of the beam 
profile), beam length a z (given in mm), and the beam energy E e . For f^j^ n (x,Q 2 ), the 
bremsstrahlung distribution function of photons in the electron, we use the Weizacker- 
Williams distribution. Finally, /*/g m (x), the beamstrahlung distribution function of pho- 



ly/e 

'e/e 



tons inside the electron, is related to f^ m (x), and is again determined by the values of 
T, a z and E e |29[| . 

As an illustration, we plot in Fig. |l] the electron bremsstrahlung distribution function 
(dashed curve) , the electron beamstrahlung distribution (dot-dashed curve) and the convo- 
lution (solid curve). For the beamstrahlung function, we take E e = 250 GeV, T = 0.1072, 
and a z = 0.12 mm, typical for a y/s = 500 GeV e + e~ LC using X-band klystron acceler- 
ating technology. We see that the distribution of electrons within the electron is strongly 
peaked at x ~ 1, and that there is a long tail extending to low fractional momentum x. 

In Fig. ||, we show the corresponding bremsstrahlung (solid) and beamstrahlung 
(dashed) distribution functions for photons inside the electron, for the same beam pa- 
rameters as in Fig. |. In this case, both distribution functions are sharply peaked at 
x ~ 0, which indicates the presence of an intense cloud of soft photons accompanying the 
beam of electrons in a linear collider. The bremsstrahlung photon distribution function re- 
mains significant out to large x ~ 1, while the beamstrahlung photon distribution function 
(which is anti-correlated with the corresponding electron distribution) completely dies off 
for large x. 

The tree level subprocess cross section for e + e~ — > ff via s-channel 7 and Z exchange 
can be found in many texts (see for instance Ref. |^0|). The 77 — > ff subprocess cross sec- 
tion can be easily constructed by crossing, along with appropriate replacements of electric 
charge and color factors from the e + e~ — ► 77 cross section given, for instance, in Ref. [j3l"f . 
In Fig. H], we show the distribution of tau pair invariant mass for e + e~ — ► t + t~ events 
at a y/s = 500 GeV LC, using the same beamstrahlung parameters as in Fig. |l|. We have 
required that pr(r) > 5 GeV and \rj T \ < 2.5. The solid histogram represents the contribu- 
tion from the a(e + e~ — * t + t _ ) cross section. The peak at m(r + r _ ) = 500 GeV originates 
in the peak of the electron PDF f e / e (x,Q 2 ) at x = 1. In addition, the peak from the 
radiative return to the Z resonance is clearly visible at at m(r + T~) = Mz- We note that 
the enhancement that would be expected to occur at x ~ due to the photon exchange in 
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Figure 1: Electron parton distribution functions in the electron from bremsstrahlung and from 
beamstrahlung, along with their convolution, for Et, eam — 250 GeV, T = 0.1072 and a z = 0.12 mm. 

the subprocess amplitude is eliminated by the px and |ry| cuts. We also plot as the dashed 
histogram the contribution to tau pair production from the 77 — * t + t~ subprocess cross 
section. It can be seen that this contribution completely dominates the tau pair production 
cross section for invariant masses below about vi ~ 250 GeV. In particular, it completely 
overwhelms the Z resonance peak, and has the potential to yield a formidable background 
for signal events with only relatively soft fermions as visible particles in the final state. 

The various 77 — ► // subprocess cross sections for SM fermions have been incor- 
porated into the event generator Isajet v7.70 j26fl , along with the photon bremsstrahlung 
and beamstrahlung distribution functions. Electron bremsstrahlung and beamstrahlung 
distributions had already been incorporated previously. 



3. Reach for SUSY in the stau co-annihilation region 

Our goal in this section is to evaluate the reach of an e + e~ LC for SUSY in the stau 
co-annihilation region by searching for stau pair production events. Most earlier studies 
have identified the reach of a LC for SUSY in the stau co-annihilation region by searching 



for e + e or pairs which would originate from selectron or smuon pair production[32 



33, 23]; see, however, Ref. for a study of stau signals at the LC. In the mSUGRA 
model at low values of tan /3, the fi, e\ and \x\ all tend to be nearly degenerate in mass. As 
tan (5 increases, the r Yukawa coupling increases, which helps to drive the tl and tr soft 
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SUSY breaking masses to values lower than their first and second generation counterparts. 
In addition, if tan (3 is large, mixing between tl and tr states reduces mf x even further. 
Thus, for low thq and large tan /3, it turns out that can be significantly lighter than 
the lightest selectrons or smuons. In this case, there should exist portions of mSUGRA 
model parameter space where a(e + e~ — > f^f^) is large, while a(e + e~ — > e^e^) and 
<j{e + e~ — > f^ifii) are kinematically suppressed or forbidden. 

For all our computations, we use Isajet 7.70 f2lf| which allows for the use of polarized 
beams, and also allows for convolution of subprocess cross sections with electron and photon 
PDFs to incorporate initial state bremsstrahlung as well as beamstrahlung effects. We use 
the Isajet toy detector CALSIM with calorimetry covering the regions —4 < rj < 4 with cell 
size Arj x A(p = 0.05 x 0.05. Electromagnetic energy resolution is given by AE em /E em = 
0.15/V-Eem©0.01, while hadronic resolution is given by AEh/Eh = 0.5/^/^00.02, where 
© denotes addition in quadrature, and energy is measured in GeV. Jets are identified using 
the Isajet jet finding algorithm GET JET using a fixed cone size of AR = VAr/ 2 + A</> 2 = 
0.6, modified to cluster on energy rather than transverse energy. Clusters with E > 5 GeV 
and \rj(Jet)\ < 2.5 are labeled as jets. Muons and electrons are classified as isolated if they 
have E > 5 GeV, \r)i\ < 2.5, and the visible activity within a cone of R = 0.5 about the 
lepton direction is less than maa;(£^/10 GeV, 1 GeV). Finally, "r-jets" are defined as 
jets fulfilling the above jet criteria, but in addition having just one or three charged tracks 
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Figure 3: Distribution in di-tau invariant mass from e+e~ — > t+t~ at a -^/s = 500 GeV linear 
e + e~ collider including brcmsstrahlung and bcamstrahlung effects, and the two-photon annihilation 
contribution, for Eteam = 250 GeV, T = 0.1072 and <r z = 0.12 mm. Cuts described in the text 
have been implemented. 



included within the jet cone. 2 

In this section, we assume an integrated luminosity of 100 fb _1 for both a y/s = 500 
and 1000 GeV e + e~ LC. We assume right polarized electron beams with ^(e - ) = —0.9, 
to minimize SM background from W pair production, with no cost and even a modest 
enhancement of the stau signal. 

We generate SUSY signal events in the mSUGRA model parameter space, plus all 
2^2 SM background processes as incorporated into Isajet. We then require both signal 
and background events with 

• no isolated leptons, 

• two tau jets. 

2 This tends to overestimate the chance that a QCD jet fakes a tau jet. We will see, however, that 
our cuts are efficient in removing this non-physics background, and the reach in channels involving tau is 
ultimately limited by the SM sources of taus. 
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The requirement of two jets restricts E v i s > 10 GeV, as each jet has E > 5 GeV. At this 
point, the signal in the stau co-annihilation region is dominated by tau jet pairs originating 
from stau pair production, while background is dominated by 77 — > t + t~ production. The 
background t + t~ pair comes out back-to-back in the transverse plane, and this correlation 
is largely maintained by the visible tau jets. The distribution in cos </>(jj), the dijet opening 
angle in the transverse plane, is shown in the first frame of Fig. |] for the t + t~ background 
(red histogram) and all other SM backgrounds {e.g. WW, ZZ production, etc.) (black 
histogram). The signal distribution is shown in the upper right frame, for the mSUGRA 
point m = 200 GeV, my 2 = 520 GeV, A = 0, tan/? = 30 and jj, > 0. We also take 
mt = 180 GeV. While the signal has some peak near cos <f>(jj) ~ —1, it maintains a broad 
distribution for all other cos <fi(jj) values. The 77 — > t + t~ background is, however, much 
more sharply peaked at coscj)(jj) ~ —1, which leads us to propose the cut: 

• cos 4>{jj) > —0.9 . 

After the cos </>{jj) cut, we show the distributions for background and signal events in 
missing transverse energy Et (middle two frames) and visible energy E v i s n,i e (lower two 
frames). The optimal cuts in these latter two quantities vary depending on where one is 
in parameter space. We examine the possibilities: 

• E T > 0, 5, 10, • • • , 195, 200 GeV, and 

• E visMe < 15, 20, 25, • • • , 495, 500 GeV. 

For each set of possibilities, we consider only cases with at least 10 signal events for 100 fb -1 
of integrated luminosity. We then pick the set of Et, E v i S M e cuts which maximizes the 
S/(S + B) ratio (S = signal and B = background) and at the same time yields a 5a signal 
which is our criterion for observability against SM backgrounds. 3 For the parameter space 
point shown in Fig. |], the total cross section before cuts is a(e + e~ — * f^f^) = 17.65 fb. 
The optimal cuts turn out to be Et > 5 GeV and 10 GeV < E vis ii,i e < 20 GeV. Using 
these cuts, the signal cross section is 0.96 fb, while the background level is 2.13 fb. We 
have checked that S/B ~ 0.3 — 1 are typical. 

In Fig. |, we show the mo vs. mi/2 mSUGRA model parameter plane for Aq = 0, 
tan f3 = 30 and fi > 0. The red shaded regions are excluded either because the fi becomes 
the LSP (in violation of bounds from searches for stable charged relics from the Big Bang) 
or due to lack of REWSB. The yellow shaded region is excluded by LEP2 searches for 
chargino pair production, and the region below the yellow contour is excluded by LEP2 
searches for a SM-like Higgs boson, assuming that this bound applies to h. We evaluate the 



neutralino relic density using the DarkSUSY program[35] interfaced to Isasugra. Most of 



3 We recognize that by trying many sets of cuts to optimize the signal as a function of parameters, we 
should include a trials factor when assessing the statistical significance of the signal. We have not done so. 
We remark, however, that current projections for the expected integrated luminosity of a LC are several 
hundred fb . Any signal identified by cuts optimized using the first 100 fb -1 of integrated luminosity, 
should therefore be visible in subsequent runs where the analysis is performed with these same cuts so that 
there is no trials factor for this new analysis. 
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Figure 4: Distribution of events in cos 4>(jj), fir and E V i S iu e for di-tau jet+ ft events from rf 
production via both e + e" and 77 annihilation, along with background from all other SM processes 
(left-hand side). We also show distributions from stau pair production signal events at mSUGRA 
point mo, mi/2, Aq, tan/3, sign(^i) = 200 GeV, 520 GeV, 0, 30, +1, with m t = 180 GcV for a 
y/s = 500 GeV LC (right-hand side). The relic density for this point is ^j^/i 2 = 0-09. 



the remaining parameter space is excluded by the WMAP bound on Q^h 2 . The exceptions 
occur in the blue shaded region, where 0.095 < f^^ 2 < 0.129, and where neutralinos can 
make up all the CDM of the universe, or in the green shaded region, where Q^Ji 2 < 0.095, 
and additional CDM candidate(s) must exist to saturate the WMAP value of QcDhih 2 ■ The 
thin region at low mo along the boundary of the excluded region is the stau co-annihilation 
corridor. The region at low thq and low is the bulk annihilation region, and generally 
gives rise to a Higgs boson mass in conflict with bounds from LEP2 searches. There is 
a remaining green region at low m^ 2 at the edge of the LEP2 excluded region where 
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neutralinos can efficiently annihilate through the light Higgs pole: ~ m^. Not shown 

is the HB/FP region at very large mo, since we want to focus on the stau co-annihilation 
region in this section. 




Figure 5: Reach of a linear collider for supersymmetry in the mSUGRA model for y/s — 500 and 
1000 GcV in the stau co-annihilation region, for tan/3 = 30, Ao = and /i > 0. The reach via 
selectron, smuon and chargino pair searches is denoted by the dashed contours, while the additional 
reach due to stau pair production is denoted by the solid contour. The dotted contour denotes the 
kinematic limit for stau pair production. The red region is theoretically excluded, while the yellow 
region is excluded by LEP2 measurements. Below the yellow contour, m/, < 114.4 GeV. The blue 
region is within the WMAP O^/i 2 2a limit, while the green region has h 2 below the WMAP 
2a limit. 



The dashed contours in Fig. || denote the reach of a y^s = 500 GeV or 1000 GeV e 



+ - 



LC for SUSY as calculated in Ref. [23]. For the -y/i = 1000 GeV runs, we use T = 0.29 and 



a z = 0.11 mm. The horizontal portion of the contour denotes the upper limit of parameter 
space which is exp lor able via chargino pair searches, while the rising dashed contour at 
lower mo values gives the reach in mSUGRA due to selectron and smuon searches. The 
dotted contours denote the kinematic limit for stau pair production at a 500 or 1000 GeV 
LC. The solid contour marks the boundary of the added region where there is a 5<r signal 
for di-tau jet events, assuming 100 fb _1 of integrated luminosity. From these contours, we 
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see that a linear e + e~ collider can see most of the additional parameter space which is 
accessible to stau pair searches. However, when m^ ~ mg along the border of parameter 
space, the tiny mass gap m-pj — yields very low visible energy, and the search contours 
turn over. Thus, in the tan (3 = 30 case shown, a direct search for stau pair production 
allows only a portion of the WMAP allowed region to be covered. We note here that 
we have also investigated stau pair signals consisting of an isolated lepton plus a single 
tau jet. In this case, no additional reach was obtained, and in some cases the reach is 
even diminished, due to the large backgrounds (mainly from WW events) in the i + r-jet 
channel. We also note here that we have investigated the lower tan (3 = 10 value. In this 
case, rrifj ~ m e R , and hardly any reach is gained by looking for SUSY in the tau pair 
channel, as opposed to the dielectron or dimuon channels. 

In Fig. |6|, we show the same tuq vs. m^2 plane, except we now take tan/3 = 45 and 
fi < 0. The shaded regions and contours denote the same constraints as in Fig. ||. The 
most noteworthy point is that the WMAP allowed region has greatly expanded compared 
to the tan /? = 30, \i > case shown in Fig. [j| This is because the A annihilation funnel has 
moved into the parameter space shown, and overlaps with the stau co-annihilation region. 
We can also see that in this case, the additional reach gained by the ditau-jet signal is 
considerable, and covers a rather large swath of the WMAP allowed region, for both a 500 
GeV LC as well as for a 1000 GeV LC. 

In Fig. [7|, we again show the mSUGRA model niQ vs. m^j plane, but now for tan (3 = 
52 and fi > 0. In this case, the excluded region where f\ is the LSP has greatly expanded, 
owing to the large r Yukawa coupling and large left-right mixing, both of which act to 
reduce the f\ mass relative to rrie R and mp -R . In this plot, the A annihilation funnel is not 
apparent (it would be located in the low mo excluded region) , but its effect is felt over 
much of the parameter space by adding to (via off-shell A and H annihilations) the stau 
annihilation and bulk annihilation cross sections, and thus enlarging these regions. We can 
see that in this case, the additional reach gained by the stau-pair search allows a ^/s = 500 
GeV e + e~ collider to fully access the bulk annihilation region at low mo and low m^j, 
which would otherwise not be accessible to searches for first and second generation sleptons. 
In addition, some portion of the stau co-annihilation region can also be explored by the 
di-tau signal, although the region with nearly degenerate f\ and Z\ remains inaccessible 
due to the tiny mass gap, and low visible energy release from tau slepton decays. Likewise, 
the \fs = 1000 GeV LC can explore all this and more, including a substantial chunk of the 
stau co-annihilation corridor, provided the mass gap between f\ and Z\ is large enough to 
yield a sufficient rate for observable signals after cuts. 

We note in passing that in the stau coannihilation region the reach of the LHC exceeds 



that of even the 1 TeV LC |23]. 



4. Reach in the HB/FP region for m t = 180 GeV 



In Ref. [23], chargino pair production was examined against SM backgrounds using a set 



of standard cuts first suggested in Ref. []32|. Briefly, the standard chargino pair cuts are as 
follows. Following Refs. and [33], it was required to have one isolated lepton plus two 
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m (GeV) 

Figure 6: Reach of a linear collider for supersymmetry in the mSUGRA model for ^/s — 500 and 
1000 GeV in the stau co-annihilation region, for tan/3 = 45, Aq = and fi < 0. The contours and 
shadings are as in Fig. |[ 



jets with, i) 20 GeV < E vis < yfe - 100 GeV, it) if Ejj > 200 GeV, then m(jj) < 68 GeV, 
Hi) EJp is > 25 GeV, iv) \m{iv)-M w \ > 10 GeV for a pair hypothesis, v) \ cos 6»(j')| < 0.9, 
| cos < 0.9, -QiCos9 e < 0.75 and Q e cos 6(jj) < 0.75, ot) 6» ac0 p(WTy) > 30°, for a 
pair hypothesis. The reach for 11 + 2j+ Et events from chargino pair production was 
evaluated using a left polarized beam with Pl = +0.9. These cuts worked well for almost 
all of parameter space, save for the large portion of the HB/FP region. In that region, 
the W\ — Z\ mass gap becomes so small that there is very little energy release in chargino 
pair production events, and the signal would usually fail the standard cuts. A specialized 
set of cuts was proposed in Ref. [23] to access the signal in the large mi/2 portion of the 



HB/FP region. These require: 

• 1 isolated lepton plus two jets, 

• 20 GeV< E msible < 100 GeV, 

• cos 4>(jj) > -0.6, 

• rn(£,j near ) > 5 GeV. 
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m (GeV) 

Figure 7: Reach of a linear collider for supersymmetry in the mSUGRA model for ^/s — 500 and 
1000 GeV in the stau co-annihilation region, for tan/? = 52, Aq = and fi > 0. The contours and 
shadings are as in Fig. |[ 

The cut on cos is on the dijet opening angle in the transverse plane. Backgrounds 

from 2 — ► 2 processes were evaluated with Isajet, while 77 — * ff background was evaluated 
using Pythia, which included only the bremsstrahlung portion of the photon PDF. The 
total SM background from 2 — * 2 processes plus 77 — * ff was found to be 0.97 fb for 
a y^s = 500 GeV LC. Using these cuts, it was found that the chargino pair signal could 
be seen essentially up to the kinematic limit even in the high m\/2 portion of the HB/FP 
region where there is no observable signal at the LHC. 

In this section, we re-examine the chargino pair signal, but this time including as 
well the 77 — > cc and bb backgrounds from both bremsstrahlung and beamstrahlung, as 
incorporated into Isajet 7.70. The background estimate from 2 — ► 2 processes plus 77 — * ff 
obtained for a <Js = 500 GeV LC with T = 0.1072 and a z = 0.12 mm is now 44.1 fb, after 
the above cuts are applied - an increase by a factor of ~ 40! (The background level using 
standard cuts hardly changes after including beamstrahlung photons, since that set of cuts 
is also adept at removing the 77 background.) 

The situation is illustrated in Fig. ||, where we show the distribution in cos6(j) for the 



-13- 



most energetic jet in signal and background events which pass the specialized HB/FP cuts. 
The mSUGRA model signal point is taken to be mo = 8850 GeV, mu^ = 400 GeV, Aq = 0, 
tan/3 = 30 and \x > 0, where we have adopted the value m t = 180 GeV for the top quark 
mass. It is seen from Fig. ||that the signal events have a broad distribution in cos 0(j), while 
the background distribution is sharply peaked at | cos9(j)\ ~ 1. Considerable suppression 
of the 77 — ► cc, bb background at little cost to signal can thus be gained by requiring in 
addition 

• \ cos 6{j)\ < 0.8 , 

which we apply to both jets in the £+2— jet events. After applying this cut, the 2^2 plus 
77 background level drops from 44.1 fb to 0.39 fb. The signal cross section for the parameter 
space point shown drops from 39.5 fb to 28.8 fb. We have re-evaluated the reach projections 
in the mo vs. mi/ 2 plane of Ref. []23| after including the updated background including 
beamstrahlung and the additional | cos#(j)| cut, and find that the reach projections suffer 
no visible change. 
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Figure 8: The angular distribution of the hardest jet in I + jj+ ft events from sparticle pair 
production at mSUGRA point mo, m 1 / 2 , Aq, tan/3, sign(fi) — 8850, 400, 0, 30, +1, with 
m t = 180 GeV for a 0s = 500 GeV LC. The relic density for this point is tl^h 2 = 0.09. We show 
background from bb and cc production via both e + e _ and 77 annihilation, along with background 
from all other SM processes. 

After the publication of Ref. [^3|, the D0 Collaboration has announced a re-analysis 
of old data on top quark pair production, using more sophisticated techniques that retain 
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more information about the event shapes for extracting the top quark mass in the I + jets 
channelQ. They now find a value mt = 180.1 ± 5.3 GeV, which has smaller error bars 
and a significantly larger value of mj than previous measurements. Combined with CDF 



measurements, the two experiments find a world average mt = 178.0 ± 4.3 GeV[36]. The 
value of mt is an important parameter in determining the location of the HB/FP region 
in the mSUGRA model[18|. We take this opportunity to update our mSUGRA reach 



projections assuming a value of mt = 180 GeV, whereas previous results used mt = 175 



GeV|3|. 



To illustrate the effect of the higher top quark mass, we show in Fig. ^ the mo vs. 
plane for Aq = 0, tan (3 = 30 and /i > 0. The main effect of the increased value of mt 
is to push the HB/FP region from the vicinity of mo ~ 3 — 8 TeV (for mt = 175 GeV) 
to mo ~ 8 — 14 TeV (for m t = 180 GeV). Another less noticable effect is that the larger 
mt value increases the relative value of rrih via radiative corrections, so that the region of 
good relic density from neutralino annihilation via s-channel h exchange has been pulled 
to larger m 1 / 2 values. This light higgs annihilation corridor is shown as a discontinuous 
narrow band around m\ii ~ 150 GeV in the figure. 

We also show the reach of the Tevatron via the isolated trilepton channel (following 
procedures listed in Ref. |[8]]), and the reach of the CERN LHC (following procedures 
listed in Ref. p2|). The LC reach is plotted as in Ref. [23]. The Tevatron reach contour 



here is plotted assuming a 5<r signal with 10 fb _1 of integrated luminosity. It remains 
qualitatively similar to the tan/3 = 30, mt = 175 GeV case, except that it is effectively 
stretched out in mo- The light Higgs annihilation corridor, which was not apparent for 
tan/3 = 30, mt = 175 GeV, is now visible, and within the reach of the Fermilab Tevatron. 
The LHC reach contour, plotted for a 5cr signal with 100 fb _1 of integrated luminosity, 
extends to my 2 ~ 1300 GeV for low mo (corresponding to a value of m g ~ 3 TeV). As 
mo increases, the contour drops until mo ~ 3 TeV is reached, whereupon it levels out 
and stays roughly constant in m^j ~ 700 GeV (corresponding to a value m g ~ 1.8 TeV). 
What is happening is that as mo increases, squarks become increasingly heavy, while m g 
remains roughly constant for a given value of m-i/2- Thus, at mo ~ 3 TeV, the squark 
contribution to LHC signals has essentially decoupled, and almost all the SUSY signal 
originates from gluino pair production, followed by gluino cascade decays. The LC reach 
plots are qualitatively similar to those presented in Ref. |23| , except that they are stretched 
out in mo until the HB/FP region is reached. As expected, this reach is mainly governed by 
the chargino mass and the mass gap between the chargino and the LSP. We note here that 
the shift in the HB/FP region to very large values of mo for mt = 180 GeV can be viewed 
favorably for SUSY theories, in that the large scalar masses will give further suppression 



to possible flavor changing or CP violating processes |37], 13]. 
5. Conclusions 

Significant portions of the mSUGRA parameter space where the neutralino relic density 
is in accord with WMAP analyses occur around the boundaries of parameter space in 
the mSUGRA region. These boundaries occur where m^ — — ► in the stau co- 
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mSugra with tanft = 30, A Q = 0, \x > 0, m { = 1 80 GeV 




m Q (GeV) 

Figure 9: Reach of the Fermilab Tevatron, CERN LHC and a linear collider for supersymmetry in 
the mSUGRA model for y/s = 500 and 1000 GeV for tan/3 = 30, A = 0, fi > and m t = 180 GeV. 
The red region is theoretically excluded, while the yellow region is excluded by LEP2 measurements. 
Below the yellow contour, my l < 114.4 GeV. The blue region is within the WMAP fig /i 2 2a limit, 
while the green region has fig /i 2 below the WMAP 2a limit. 

annihilation region or where tu^ — — > in the HB/FP region. In these regions, pair 
production of NLSPs will result in collider events with very low energy release, and low 
invariant mass. The process 77 — * ff where the 7s come from initial state bremsstrahlung 
or beamstrahlung is also highly peaked at low // invariant mass, and may thus constitute 
an important background process in these DM allowed regions of parameter space. 

We have incorporated 77 — * ff into the Isajet event generator for e + e~ collisions. 
This enables us to evaluate signal and background rates for e + e~ — ► f^f^ — > tt+ fix 
production in the stau co-annihilation region. Using suitable cuts to tame the 77 — * tt 
background, we find that significant regions of additional reach are obtained, especially if 
tan j3 is large. 

We have also investigated the 77 — ► cc, bb background to e + e~ — ► W^W^~—> P+jj+JSr 
production in the HB/FP region of the mSUGRA model. Inclusion of the beamstrahlung 
backgound significantly increases the SM background for this signal. However, the addi- 
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tional background can be effectively removed by requiring the jets to be away from the 
beam directions. 

We update our mSUGRA reach plots for tan j3 = 30 by increasing the value of mt from 
175 GeV to 180 GeV, in accord with recent D0 measurements. This has the effect of 
pushing the HB/FP region out to much large values of mo, and effectively stretching out 
the reach contours for the Fermilab Tevatron, the CERN LHC and an e + e~ LC. Expressed 
in terms of the chargino mass, the reach in the HB/FP region is qualititatively unaltered. 
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